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SUMMARY 


By the vise of methods of linearized nonviscous, conical— flow 
theory, formulas for computing the rolling and yawing moments for flat, 
thin, sweptr-hack wings at angle of yaw have Been developed. The for- 
mulas were derived specifically to treat that family of wing plan forms 
for which all edges are subsonic (i.e., the component of the free— stream 
velocity normal to the edge is subsonic) and straight. The formulas are 
also applicable to wings for which all but the trailing edge is sub- 
sonic. The rolling and yawing moments of several representative plan 
forms were calculated. All the wings showed stable rolling moments 
about the axis of symmetry of the wing plan form, and showed stable 
yawing moments about a vertical axis through the apex of the wing. 


INTRODUCTION 


In previous treatments of the problem of determining the rolling 
or yawing moment due to steady sideslip or yaw at supersonic speeds 
(references 1, 2, 3* and. 4), the case of thin, flat, swept— back wings 
with all edges subsonic (fig. 1) — that is, wher$ the component of the 
free-stream velocity normal to any edge is subsonic - has not been 
treated. The classes of swept-back wings studied in the references 
cited have supersonic trailing edges, where the Kutta condition is 
fulfilled without affecting the existing flow over the wing. However, 
in fulfilling the Kutta condition at a subsonic trailing edge the flow 
over the wing will be affected. The Kutta condition may be realized 
in such cases if the lift distribution (i.e., streamwise component of 
the perturbation velocity) is made to go continuously to zero at these 
edges. Lagerstrom in reference 5 suggested that this may be accomplished 
by the superposition of conical-flow fields on semi -inf inite flow fields 
to cancel lift distribution beyond the edges. The methods suggested in 
reference 5 were used in reference 6 to calculate the lift and pitching 
moments on swept-back wings having all edges subsonic. In a similar 
manner the methods of reference 5 are utilized in the present paper to 
obtain the rolling moments of such wings at an angle of yaw. 
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In addition to satisfying the Kutta condition at the trailing edge, 
it must also he satisfied for the trailing tip in the case of the yawed 
wing. In figure 2, it can he seen that as the wing is yawed, one initi- 
ally streamwise tip becomes a- subsonic trailing edge where the Kutta 
condition must apply, the other tip becomes a subsonic leading edge where 
the lift distribution is similar to that associated with the leading edge 
of a thin flat plate in subsonic flow. At such leading edges the linear- 
ized theory gives an infinite value for the streamwise perturbation 
velocity (i.e., lift distribution) accompanied by a horizontal suction 
force, These conditions at the tips are taken into account in the equa- 
tions that are subsequently derived for the lift, rolling moment, and 
yawing moment . 

Although the formulas in this paper were derived for wings with 
initially unraked tips, they may be applied to wings with initially 
raked tips if a change of the parameters involving the slopes of the 
tips is made. It should be pointed out that the rolling up of the 
vortices around a streamwise or nearly streamwise edge is not predicated 
by the linear theory. Therefore, the results given in this paper may 
not completely describe the nature of the flow at any edge which is 
nearly alined with the free— stream direction. 

SYMBOLS 


V 

M 

P 

P 

q 

Ap 

qa 

a 

S 


General 

free— stream velocity 
free— stream Mach number 
VM2-1 

free— stream density 
free— stream dynamic pressure 
coefficient of local lift per unit angle of attack 

angle of attack, radians 

angle of yaw, degrees (positive as shown in fig. 3) or negative 
sideslip 

3 tan Jr 
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A 


co 

c t 

s 

S 

A 

\ 


x,y,z 


x a^a 


x -b^y-b 


1 


Wing Parameters 
/ 


aspect ratio 

• \ % 

root ohord 

tip chord 

semi span 

area of wing plan form 

angle of sweep of leading edge, degrees 

taper ratio 

Coordinates 

/ 

Cartesian coordinates in the stream direction, across the stream., 
and in the vertical direction, respectively 

coordinates of the apex of a constant— load element defined by 
ray a 

coordinates of the apex of an element used in secondary correc- 
tions 

coordinates of the apex of the trailing edge 

spanwise coordinate referred to the apex of a constant-load 
element 



In the following, all slopes are measured in the positive sense 
counterclockwise from the x axis: 


a 


t a 


% 


*5 


P x (slope of ray from apex of wing) = p — 
subscript) 


P x (slope of ray through x a ,y a ) = p 


p x (slope of ray through x- b ,y b ) = p 


y-y a 

x-x a 

y-yb 

X-X-u 


P x (slope of ray through Xg,yg) = p 


y-y& 

X-X c 


(used also as a 




k 
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Rays Describing Wing Geometry 


The following are for the unyawed wing shown in figure 1: 
a t p x (slope of ray from apex through intersection of right-span 

tip and trailing edge) .^ s > . 

c 0 +(P&/®t) 

m p x (slope of right-span leading edge) = p cot A 

mj. p X (slope of right-span trailing edge) 


The following are for the yawed wing shown in figure 2. Here the 
subscripts l and 2 -refer to the right and left spans, respectively. 

a i. > a 2 + 3 x (absolute value of the slope of the ray from apex through 
the intersection of the tip and trailing edge) 



mi,m 2 P x (absolute value of the slope of leading edge) 

P X (absolute value of the slope of the trailing edge) 



In general, any yawed rays di, d 2 which describe the wing geom- 
etry-may be obtained from the unyawed rays d by the following rela— . 
tions : 

Perturbation Velocity Components 

*, ' 

u,v,w perturbation velocities in the x,y,z directions, respectively 

u^(a) streamwise perturbation velocity for a triangular wing at angle 

of yaw , 
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u Aj(a) value of %(a) for wing reflected in the x axis 
Ua incremental velocity on any constant— load element 

u value of Ua for wing reflected in x axis 

Ujj incremental velocity on the elements used in secondary correc- 

tions 


Uc 


value of u« (a) for a = — 5 


Forces and Moments 


DL ' 

Lo 

AL 

(AL) a 

V 

AL’ 

(AL')a 

N 

a 

c d l 

C L 

Cl 


'“o 


drag due to lift 

"basic lift for the entire wing 

correction to "basic lift 

correction to basic lift due to the application of one constant- 
load element 

"basic rolling moment for the entire wing evaluated about the 
axis of symmetry 

correction to the basic rolling moment 

correction to the basic rolling moment due to the application of 
one constant— load element 

yawing moment about the apex of the wing 

suction force per unit length of edge 

drag coefficient^ 


lift coefficient 


tt?) 


rolling-moment coefficient ^ rolling _ moment ^ 

f yawing moment 


yawing-moment coefficient ^ 2qsS 

O 

coefficient of yawing moment due to normal force (ocCj) 
coefficient of yawing moment due to suction force 
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Elliptic Integrals 


k modulus of elliptic integral (value given in text) 

k' complementary modulus ( J 1— k 2 ) (value given in text) 

K(k) 1 complete elliptic integrals of the first and second kinds, 
E(k ) J respectively, of modulus k 

F(cp,k) \ incomplete elliptic integrals of the first and second kinds, 
E(cp,k)J respectively, of modulus k and argument cp 

A 0 (tp) - [E(k)F(k',cp) + K(k)E(k,q>) - K(k)F(k' ,q>) ] 

K 

ANALYSIS 


The present report utilizes the same cancellation— of— lift process 
as was used in reference 6 to determine the induced effects due to the 
interaction between the flows on the upper and lower surfaces of wings 
with subsonic tips and trailing edges. The method consists of super- 
posing conical— flow fields on a semi— inf inite flow field in order to 
satisfy the boundary conditions about a given wing plan form. The analy- 
sis of the characteristics of a thin, flat, yawed, swept wing with sub- 
sonic edges by this method falls logically into two main divisions: 
namely, (l) that of computing the basic forces and moments on a given 
plan form due to the flow field associated with a semi— inf inite yawed 
lifting triangular wing of the same apex angle as the given plan form, 
and ( 2 ) that of determining the forces and moments induced by the inter- 
action between the upper— and lower— surface flows around the subsonic 
tips and trailing edges. These induced forces and moments oh the wing 
may be computed by superposing conical— flow fields beyond the tips and 
trailing edges in order to cancel the lift of the. infinite field. 


The flow fields which are superposed in the cancellation process 
may be expressed in terms of the perturbation velocities (u,v,w) . If 
Q represents a perturbation velocity it must satisfy the linearized 
equation for, a nonviscous, irrotational compressible fluid, namely. 


(1-M 2 ) ~ + 


cfQ 

cbr 


d 2 Q S_Q 
dy 2 dz 2 


= 0 


CD 


It can be shown that Q may be expressed as a function of a single 
complex variable; the variable is homogeneous of degree zero, that is, 
constant along rays from a fixed point in space (references 5 and 7). 

The flow fields represented by Q therefore are conical. : Several of 
these flow fields are described in detail in reference 5. The boundary 
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conditions which will be imposed on these fields for purposes of the 
present problem will be indicated later. 

In order to simplify the analysis, the wing plan form has been 
divided into a number of separate regions as shown in figure 2. The 
regions are defined as follows: 

1. The basic region, is the entire wing plan form. In this region 
first approximations to the lift and rolling moment about the axis of 
symmetry are obtained by integrating over the entire plan form the lift 
distribution given by the expression for a yawed triangular wing having 
the same apex angle as the given plan form. This lift and rolling 
moment will be designated as the basic lift and basic rolling moment. 

2. The primary regions are the regions which are influenced by the 
cancellation of the basic lift distribution for a sideslipping triangular 
wing beyond the tip and trailing edges of the desired plan form. The 
corrections to the lift and rolling moment which occur in these regions 
will be called primary corrections . 

3. The secondary regions are those which are influenced by addi- 
tional corrections which must be applied at the edges of the wing where 
the cancellation flows introduce residual lift distributions beyond the 
confines of the wing. These corrections to the lift and rolling moment 
will be called secondary corrections, and, in general, are negligible 
for small angles of yaw. 


Basic Lift and Rolling Moment 


The basic lift is obtained by integrating the lift distribution of 
a yawed semi— infinite triangular wing over the wing plan form shown in 
figure 2. If mx/p and m 2 /p sure the magnitudes of the slopes of the 
right- and -left— span leading edges of a yawed triangular wing with sub- 
sonic lea.ding edges, the velocity distribution in the stream direction 
may be obtained from reference 8 as 


P 0 Va[2m 1 m 2 +a(m 1 -m 2 )] 

u A (a) = . (2) 

Pv^ (mx-a) (m 2 +a) 


where 


Po = 


1 

E(ko) [ V(l+mx)(l+m 2 ) + J ( 1-mi ) ( l-m 2 ) ] 


(3) 
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and E(k 0 ) is the complete elliptic integral of the second kind of 
modulus 


*o = 


2[ ( 1— ffll 2 ) ( ^ 

4 ( 1+mi) ( l+m 2 ) + J ( 1-mi) ( l-m 2 ) 


W 


The velocity distribution u.(a)' at a wing section normal to the flow 
direction is shown in figure j for typical values of m and ,8. The- 
corresponding lift distribution per unit angle of attack for the tri- 
angular wing may be expressed as 


AP = ^(a) 
qa Va 


(5) 


Integrating equation (5) over the wing plan form shown in figure 2 
yields for the total basic lift (see appendix A) 


P 2 +8 2 


^ = 2P 0 s 2 „ 
qa P 2 


■ Il(B) ^^) 2 + I 2 (8) + (^§^ I 3 (5KI 4 ( 5 ) 


a 2t -6 J 


. ( 6 ) 


The values of Ij.(&) 'i=l,2,3A are given in appendix A (equations 
(A8) and (A9)). ' 


The basic rolling moment about the axis of symmetry of the wing is 
given by 



+ 


I 2 (8) - 


113 ( 5 )^ 


m 2 t -a 2t \ 3 

a 2 t — 5 J 


I 4 (&) 


( 7 ) 


See appendix A, equation (A12), for IIi(5) and 113 ( 8 ). 
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Primary Corrections 


Tip corrections — right span (region II) .— The following correc- 
tion applies over the region shown in figure 2 as ABF, where AB 
corresponds to the right-span tip. This tip acts like a subsonic trail- 
ing edge since the streamwise velocity component must cross this edge to 
leave the wing. It is therefore necessary to impose the Kutta condition 
along ^this edge. This condition is satisfied by canceling the lift dis- 
tribution due to the triangular wing beyond the confines of the tip 
edge. The cancellation is accomplished by superposing the flow fields 
of an infinite number of constantly loaded overlapping sectors of infin- 
ite extent along the tip. As shown in figure 4, each of these sectors 
is bounded by a line along the wing tip and a fixed ray a = constant 
from the apex of the wing. It is prescribed that each sector carry a 
constant load proportional to the incremental change of the loading for 
the triangular wing along the rays, namely. 


du A (a) 
u a = — ^ 


da 


( 8 ). 


The effect on the wing of one constant— load sector, defined by a 
single ray a, is first computed, and then the integrated effect of all 
the sectors along the tip is determined. The solutions for constant- 
load elements should be such that they introduce no change in the angle 
of attack of the wing j that is, the induced downwash on the wing should 
be zero. In the stream the pressure difference must vanish. The 
boundary conditions for a given elemental sector may be expressed as 
follows (fig. 4): 


1. u = Ua — 5 5; t a < a 

2. u = 0 a < t a — + 1 

3. w = ,0 -1 < t a < - 5 

An additional condition on the perturbation velocities is that they 
must vanish on the Mach cone from x a , y a . If a general ray from the 
apex of one of these constant-load elements is defined by 


ta - 0 


y^y a 

x— x a 


(9) 


where 
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a = P ^ 
x n 


( 10 ) 


then the function satisfying the conditions specified, given hy 
Lagerstrom in reference 5, may he expressed as 


i 


u = r.p 


tli(a,t a ) = 


* a® C0S ~ 1 

(a+6)(l+t a )+(t a +8)(l+a) 

(t a -a) (l-r5) 


(ID 

( 12 ) 


The change in local pressure for a single constant— load sector may 
therefore he expressed as 


Ap\ _ b Ua 
qaj ~ it Va 


cos -1 Tji(a,t a ) 


r (13) 


-1 < t a < - 5 


The integration of the lift for a single element in region II gives 


(AL) a = pVua s 2 Gx(a) 

P 


( 1 *) 


where Gi(a) is given in appendix B (equation (Bll)). 

The total induced lift correction on this tip is determined hy 
calculating the effect of canceling with one sector the lift density at 
the leading edge and then adding hack the integrated effect of all the 
elements along the tip edge. Since the magnitude of the lift density at 
the leading edge is infinite, the quadrature described above is performed 
as follows: 


AL = pVs 2 


S 2 +5‘ : 


lim 

a — •> mx 



a)Gx(a) 



du A (a) 


Gx(a) da 


3 


da 


(15) 
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Integrating by parts and passing to the limit gives 



Similarly, the rolling-moment correction for a single element is 


2 e2 

(AL’ ) a = - PV u a s 3 £-jp- H^a) (17) 


and the total induced moment may be expressed as 


AL 1 Vs 3 ( ft 2 +S g ) 

qa qa(3 


The derivatives Gi'(a) and Hi’ (a) 

(B13) and ( B22 ) ) . 

. Tip correction - left span (region III) .- The following correction 
applies over the region shown in figure 2 as DGC and includes the left 
tip DC which may be considered a subsonic leading edge, since the 
streamwise component of the velocity must cross this edge in order to 
reach the wing. 

In canceling the lift distribution beyond a leading tip, it is 
necessary to modify the previously used cancellation function (equation 
(ll)) to satisfy the condition of uniform downwash over the wing plan 
form (reference 5). The modification takes the form of an additional 
term, which in turn introduces a singularity in the lift distribution at 
the leading tip analogous to that at a leading edge in subsonic flow. 

In treating the left-span tip, it is convenient, as in appendix A, to 
consider the wing to be reflected in the x axis so that 'the left-span 
tip may be treated in the first quadrant. In figure 5 the boundary con- 
ditions for a single element and the regions of influence for the primary 
induced corrections are given for the reflected wing. The boundary con- 
ditions are as follows: 

1. u = Ug^ 5 < ,t a < a 

2. u = 0 a < t a < 1 


r i 

u A ( a )Hi’ (a) da 

a l. 


( 18 ) 


are given in appendix B (equations 


3- 


w = 0 


- 1 < t a < 5 
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The value of u^ is the value of Uq, for the wing reflected in the 
x axis. (See appendix A, equation (Al), for the required transforma- 
tion.) 

The function which satisfies the boundary conditions along a lead- 
ing tip is givep in references 2 and 5 and takes the form 



25 

a( 1+5) 



(a— 5) ( 1-fa) ( 1+ta) 
5-t a 


- 1 < t a < 5 


where 


Tl2(a,ta) 


(a— 5) (l+t a )+(t a — 5) (l+a) 
("t a — a) (1+5) 


(19) 


( 20 ) 


The lift distribution for one element becomes 


d 



ku *l 

itVa 


cos -1 Ti 2 (a,t a ) 


25 / ( a— 5 ) ( l+a) ( l+t a ) 

a(l+5) v S— ' t a 


( 21 ) 


The integration of lift for a single sector gives 


32+52 

(AL) a = pVu a s 2 - G 2 (a) (22) 

1 P 

The total lift correction is then 
AL B 2 +5 2 

— = - pv s2 — — J u Az (a)G 2 ’ (a) da (23) 

a 2t 

where the functions G 2 (a) and its derivative G 2 ' (a) are given in 
appendix B (equations (Bll) and (B13), respectively). 

The rolling moment induced by a single sector becomes 


(AL')a = pVu ai s 3 H 2 (a) 


(24) 
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and the total rolling-moment correction takes the form 


fr - - p sa P " A ! (a) H =' (a> 48 (25) 

a2 t 

The function H 2 (a) and its derivative H 2 ' (a) are given in 
appendix B (equations (B20) and (B22), respectively). 

Wake corrections (regions IV and V) .— The following corrections 
apply over the regions defined in figure 2 as O'EAB and O' Cl, where 
O'B and O'C are Both subsonic trailing edges at which the Kutta condi- 
tion must be imposed. Behind these edges, the lift distribution is 
reduced to zero in two steps (fig. 6): (l) canceling lift in the amount 
of the lift density at the apex of the trailing edge with a single semi- 
infinite field extending over the whole wake region, and (2) removing 
the remaining variation by the superposition of infinitesimally loaded 
elements along the trailing edges in the manner described for the tip 
edges. The latter elements are bounded by the trailing edge an d rays 
from the apex of the wing. 

The correction to the lift resulting from step.(l) above for the 
single field influences both regions IV and V. If 


t& =• P —S (26) 

x x 8 

refers to a ray from the apex Xg, y 5 of the trailing edge, the boundary 

conditions which must be satisfied for the correction function are 

1. u = u a (-5) - m 2t < t 5 < mi t 

2. w = 0 - m 2t > tg > m lt 

These boundary conditions and the regions affected by the corrections 
are shown in figure 7. The function takes the form 

u(t 5 ) = r.p. ug (27) 


where 
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The foregoing expression is a generalization of the equation given in 
reference 9 for the symmetric trailing-edge correction and is valid when 
the wing is yawed. The u§ is that given "by the triangular-wing solu- 
tion for a equal to — 5, or 


u 5 = V- 6 ) = 



(30) 


The total correction to the lift induced by the cancellation of the 
field of uniform lift distribution in the wake is obtained by integrat- 
ing equation (27) over the regions of the wing contained within the Mach 
cone from the apex of the trailing edge. If the Mach lines intersect 
the leading edge, as shown in figure 2, it is convenient to integrate 
first over the entire region O' BA’ between the trailing edge and the 
Mach cone and then to subtract from this the integral over the region 
EAA’ between the leading edge and the Mach cone. This procedure is 
indicated in equation (B24) . Carrying out the prescribed operations 
gives a result in the form 


AL 

qa 



(31) 


The values of R 0 i, Rii, and R 2 i are given in appendix B (equations 
(B29 ) } (B32), and (B33))« The index i=l,2 corresponds to the right 
and left span, respectively. The terms Rqi and R02 constitute the 
lift corrections for the region bounded by the Mach cone from the apex 
of the trailing edge, the trailing edge, and the lines which define the 
tip edges. The additional terms Rii and R 2 ^ are those required if the 
Mach lines from the apex of the trailing edge cross the leading edge. 

As the angle of yaw is increased, the effect of the terms R X1 and R 21 
increases while that of the terms R X2 &nd R22 decreases. 



NACA TN 2262 


15 


The correction for the rolling moment about the axis 
resulting from the lift correction of equation ( 27 ) is 


of symmetry 


AL' 

qa 


I '->>■ {-»■.* 

i=l 


fo! 5 _ [» 1 -(-l) 1 5i 2 

S 3 (p 2 + 62) 3 P 21 S 3 (p 2 + 62) 3 


P 4 



(32) 


The value for R ai is given in appendix B (equation (B4l) ) . As in the 
lift correction, the terms Rii, R 2 i, R 3 i apply only if the Mach cone 
from the apex of the ving crosses the leading edge. 

It may be seen in figure 6 that the correction effected by super- 
posing one constant— load sector in the wake cancels the largest part of 
the triangular-wing lifting pressure behind the trailing edge. It is 
necessary to remove the remaining lift by superposing additional 
constant-load sectors at the trailing edge. The amount of lift to be 
canceled by these elements is shown in figure 6 where it is noted that, 
along the right span, some of the elements will add to the lift distri- 
bution where the amount canceled by the single field is excessive. One 
of the elemental sectors, which will be referred to as "oblique" sectors 
is also shown in figure 6. 

Considering first the right-span trailing edge, it is found that 
the function necessary to cancel the lift distribution downstream of 
this edge must satisfy the following boundary conditions for each sector: 


1. 

PS 

II 

& 

a — ^a ^ m it 

2. 

o 

II 

2 

— 1 < t a < a 

3- 

o 

II 

> 

< t a < + 


These boundary conditions and the region of the wing affected by this 
correction are shown in figure 8. 

Note that while the condition of zero induced downwash may be 
specified for the wing area adjoining the element, it is not possible to 
satisfy this condition on the opposite wing panel if the Mach cone trace 
from the apex of the element should intersect the opposite trailing edge. 
However, it has been found for the wings investigated in this paper that 
the downwash induced by the oblique trailing— edge corrections was negli- 
gible. 
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The function satisfying the specified boundary conditions for the 
oblique trailing— edge corrections is 


where 


u 


r.p. 


u A (a) 


COS - * 


H^a^a) 


( 33 ) 


B 3 (a,t a ) 


( 1-a) ( t a -m lt )-(m lt -a) ( l-t a ) 
(l-ffl lt )(t a -a) 


mi t £ t a £ 1 




(34) 


The induced lift for a single element then becomes 


ft 2 +5 2 

(AL) a = pVu a s 2 — — G 3 (a) 
P 


( 35 ) 


and the total induced lift is 


i -- » ^ Mi) a. 

-—ft 


( 36 ) 


where. G 3 (a) is a function given in appendix B (equation (Bll)). 

The rolling-moment correction due to a single constant— load sector 
is 


(AL^a = f Vu a s3 JLj±_ H3( a ) (37) 

j p 


and the total rolling-moment correction becomes 


Mil H. (a, da 

3 qap - J da 3V ' 

—5 

(For % (a), see appendix B, equation (B20).) 


( 38 ) 
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The corrections for the lift distribution, roiling moment, and lift 
for the left-span trailing edge may be obtained from the solutions for 
region IV by reflecting the wing in the x axis, that is, replacing 5, 
m i.£> and % ms t’ res P ec tively. The functions Ga(a) 

and H 3 (a) are replaced by G 4 (a) and H 4 (a) . The latter two functions 
are given in appendix B (equations (Bll) and (B20)). 


Secondary Corrections 

In canceling the lift about a given swept>-back-wing plan form in 
order to establish the correct boundary conditions, the primary correc- 
tions previously described introduce errors that can influence the lift 
and rolling moment j however, only in the case where the Mach lines from 
the apex of the trailing edge cross the leading edge are these secondary 
corrections of any considerable magnitude. In this case, the primary 
correction which may cause sizable errors is the initial wake correction 
introduced by the single field from the trailing-edge apex. This wake 
correction adds negative lift outboard of the tip station and ahead of 
the leading edge. The secondary regions affected by canceling this 
negative lift are shown in figure 2 as EABE' and II* C. The method for 
making these secondary corrections to the lift and rolling moment is 
- given in appendix C. 


YAWING MOMENT DUE TO SIDESLIP 


In order to compute the yawing moment for a thin, flat wing at 
angle of incidence in a frictionless flow, it is necessary to find the 
moment of the total force (x— y plane) about a suitable vertical axis. 
For a swept wing with subsonic leading edges, the total force in the 
-plane of the wing is the resultant of the drag force due to lift and the 
thrust force due to the edge suction. The drag due to lift is <xL for 
thin, flat wings, and, in coefficient form, may be written as 


Cd l = ojCl = 


acp 

da 


a2 


( 39 ) 


where a is measured in the plane of symmetry of the wing. The rolling 
moment about the axis of symmetry may be expressed in terms of the yaw- 
ing moment due to the drag force about an axis no rmal to the axis of 
symmetry as follows: 


L' 


T — aL _ Dl — 
= L y = — v = — u v 


a 


a 



(itO) 
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where _N 0 is the yawing moment due to the drag force Dj,. In this equa- 
tion y is the moment arm in. the plane of the wing from the center of 
pressure to the axis of symmetry of the wing. In coefficient form, the 
yawing moment due to the drag force may be written as 



a q 2 


(41) 


The suction forces are those which are associated with the singu- 
larity in the lift distribution which occurs at subsonic leading edges 
of thin, flat wings. In reference ID, the formula for finding these 
suction forces is given for any subsonic leading edge. By considering 
the flow two-dimensional in the neighborhood of the leading edge, it is 
possible to determine the formula for the suction force (per unit length) 
normal to this edge as 




(42) 


where f(x) is the strength of the leading— edge singularity in Ua (a) 
and m is 3 times the local slope of the leading edge. 

For the yawed triangular wing 


fi(x) = ( p2+ml2) 

P 


1/4 


Va 


«/ n l+ m 2 


E(*o) V ( 1+Wl ) ( l+m2 j" + (l - ®! ) ( l - ®2 ) 


Jx (43) 


where 1=1,2 refer to the right and left spans, respectively, and k 0 
is as given in equation (4). This value for ff(x) above is the same 
for the swept wing unless, the Mach lines from the apex of the trailing 
edge intersect the leading edge. For small angles of yaw the effect 
introduced by the Mach lines crossing the leading edge is negligible and 
not considered here. The yawing moment about the apex due to the suction 
on one edge is 


= irp Vl-nij 2 / ^ *fti g J' i 


fi(x) 2 xdx 


( 44 ) 


where x Qi is the chordwise distance to the wing tip from the apex. 
The total yawing moment due to suction for steady sideslip or yaw is 
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N cr - N — N <j 2 


( 45 ) 


The integration of equation (45) yields in coefficient form. 




a 2 6E(k 0 ) 2 8 3 L l+mims + y ( 1-mi 2 ) ( l-m 2 2 ) ] '“ l2) X 

(^) - ^ <^ 2) (^) 3 ] <*> 


The total yawing moment about the vertical axis (z axis) through the 
apex of a swept or triangular wing with subsonic edges is 


c n _ Ciip-Cng 
a 2 a 2 


(47) 


It will be noted that the suction effect along the tip which becomes a 
leading edge has not been treated here. The contribution of this effect 
to the ove3>-all yawing moment is small by comparison with that of the 
subsonic leading edges, and there is some question as to whether the tip 
suction is realized in a real fluid. 


APPLICATION AND DISCUSSION 


The analysis in the foregoing section has been directed mainly 
toward finding the rolling moment for a swept— back wing with all edges 
subsonic. The representative wing chosen for analysis (shown in fig. l) 
is swept back 63°, has a taper ratio of 0.25, and an aspect ratio of 4. 

In addition, the rolling moments for three other plan forms are computed 
for the sake of comparison. These wings have the same sweepback and 
wing area as the tapered swept— back wing above and are of relatively low 
aspect ratio. These plan forms are shown in figure 9 and represent a 
triangular or delta wing, a delta wing with streamwise tips, and a so- 
called "arrow wing" with a supersonic trailing edge. In figure 10, the 
variation of rolling-moment coefficient with angle of sideslip is plotted 
for the various plan forms considered. All the wings gave a positive 
rolling moment for a positive angle of yaw, that is, the left tip tended 
to move upward and the right tip downward. 



20 


NACA TN 2262 


In determining the rolling moment of the 63° swept wing, it was 
found necessary to consider only certain of the primary corrections to 
the "basic rolling moment. The primary corrections which contributed 
noticeably to the over-all rolling moment include the tip corrections 
and the wake correction due to the single uniformly loaded element. The 
magnitudes of these corrections as compared with the total rolling 
moment per degree of yaw are given in table I . The secondary corrections 
considered were found to be negligible except for the correction to the 
primary wake correction at the leading and tip edges. This correction 
is listed in table I as secondary effects. Ih the case of the delta 
wing with streamwise tips, the only corrections necessary were the pri- 
mary tip corrections. For the arrow and delta wing, it was necessary 
only to compute the basic .rolling moment. The rolling moments for these 
two wings are in agreement with the results obtained in reference 1. 

In general, when computing the rolling moment for a highly tapered 
swept wing by the superposition method, it is necessary to consider only 
the basic rolling moment and the primary correction given by equation 
(32). For the wing analyzed here, the error in neglecting all other 
corrections is less than 4 percent for the range of yaw angles consid- 
ered. It should be pointed out, however, that the primary and secondary 
tip corrections become increasingly important as the taper is decreased. 

The yawing-moment coefficients of the wings analyzed are given as 
a function of yaw angle in figure 11. The yawing moments due to normal 
force and leading— edge suction sire computed about the vertical axis 
through the apex of the wing. All wings exhibit stable yawing moments, 
that is, for a positive yaw angle the wing tends to yaw to the left or 
the left tip tends to move forward sind the right tip rearward. The 
yawing moments for wings 3 and 4 were in agreement with the results of 
reference 4. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Oct. 30, 1950. 
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APPENDIX A 

FORMULAS FOR THE INTEGRATED BASIC LIFT AND ROLLING MOMENT 
LIFT FOR THE BASIC WING PLAN FORM 


The total basic lift is obtained by integrating the lift distribu- 
tion given in equation (2) over the plan form defined as region I in 
figure 2. The integration is performed over four separate areas shown 
in figure 1 as * OBA, OBO ' , 0C0*, and OCD. For the left span a reflec- 
tion of the wing in the x axis permits the integration to be performed 
in the first quadrant. Such a reflection corresponds to the following 
transformation of parameters: 


— 6 — 5> 8 &1 t — ** ^ m 2 m 2 — ^ m i 


m^ — -5> m 2 ^. m^t 


Then the lift may be expressed in the following form: 


qa Va 



i=l 


a ii ' 


o 




where the Jacobian has the value 


and 


dj*,y) = B X 

d(y,a) a 2 


u Ai( & ) = u Ae( a ) = u A< a ) 
u A 3 ( & ) = %( a > = 


(Al) 


(A2) 


(A3) 


m 


The limits yi(a) are the y coordinates of the intersections of 
the rays from the apex and the tip or trailing edges; y^a) and y (a) 
are at the trailing edges while y 2 (a) and y 4 (a) are at the tip edges 
The values for the limits of integration are 
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/ 


y x (a) = c 0 a(mi t +5) [(m lt -a) 1 

y 2 (a) = as [(3(a+6)] */p 2 +5 2 


y 3 (a) = c 0 a(m 2t -6) [ (m 2t -a) 1 

y 4 (a) = as [p(a-5)] 1 */P 2+52 




>(A5) 




and, for the integration with respect to a, 


a n = — 5 ai 2 - a 2 i - ai^. 

a l 3 = 5 a 14 = a 23 = a 2t > (A 6) 


8-22 “ 


&24 == 


It can he seen that the third and fourth terms of equation ('A2) may 
he obtained directly from the first and second terms hy applying the 
transformation given in equation (Al) . 

The integration of equation (A2) gives for the total basic lift 


c ■ if[(=lgr) «» * '■»' *(=g=?) V« * <•»>] 

(A7) 


where 


1 r [m 2 (m lt -m 1 )-m 1 (m lt +m 2 )]y (m 1 -a)(m 2 +a) 

(mi^-mi) (m 1 ^.+m 2 ) \ 


m i t ~ a 


mi^(mi+ni 2 ) (mi— a) (m 2 +in.i^)— (1112+Q-) (mi^.— mi) 

slur 1 


2 J (m!^— mi ) ) 


|m lt -a| (m x +m 2 ) 


a_a it 

L- 8 


(A8) 


and 
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I Z (6) i -{+ -- ■ cosh- 

(mi+S) (m 2 — 8) L 2 ^(mi+b) (m.2— 6) 


—i (mi— a) (m2r-8)+(m 2 +a) (mi+5) 
(a+ 6 ) (mi+m 2 ) 




a+5 J a=a lt 


(A9) 


The values for I 3 (S) and 14 ( 8 ) may be obtained by applying transfor- 
mation (Al) to equations (A 8 ) and (A9)> respectively. 


INTEGRATED BASIC ROLLING MOMENT 


The basic rolling moment is obtained by integrating the moment due 
to the basic lift distribution about the axis of symmetry of the wing 
over region I of figure 2. The integration is performed over the four 
areas of region I as in the case of the basic lift. The basic rolling 
moment, therefore, is defined as 


hiL 

qa 


Va */p 2 +S 2 


i r *•> f 


i=l a ii 


<Ky,a) 


i r 1 <-> - - 

i =3 a ii 0 


(A10) 


Here the Jacobian haB the. same value as that given in equation (A3) and 
the limits of integration are as given in equations (A5) and (A 6 ) . In - 
the same way as in the lift case, the third and fourth terms of equa- 
tion (AlO) may be obtained from the first and second by equation (Al) . 
The integration of equation (AlO) yields for the rolling moment 



(All) 
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Here I 2 (,&) and I 4 (&) may be obtained from equation (A9) while the 
value for II 1 (b) is 


IIi(&) = - Ii( 8 ) - 


m lt +& (• 

/ m x 

m 2 \ 

‘ 1 

2 1 


m l t +m 2 J | 

_ mi t~ a 

1 


2(m 2 -mi) 

.1 i 


2 \ n ll^ +m 2 


(mi^nii) (mi^.+m 2 ) J m^-a 


*J_ (mi— a) (m 2 +a) 
t“ 


mi 


h +m 2 mi t -mi 


) 


{(ms-mi) (— — + — — (-*- 
L \mi t +m 2 mi t -^ai/ 2 \m lt ' 

[ i ♦ I ♦ *£2fV| 1 x 

2 \m lt +m 2 m^-nii/JJ 

1 , (mi.+m 2 )(m 1 -a)-(mi.-m 1 )(m 2 +a) 

sin -1 ' ' * 


J (m lt -mi) (mi t +m 2 ) 


|mi t -a| (m!+m 2 ) 


a= ai t 

a=— 5 
(A 12 ) 


The value for II 3 may be obtained from equation (A12) by equation (Al) 
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APPENDIX B 

FORMULAS FOR THE PRIMARY INDUCED LIFT AND 
ROLLING-MOMENT CORRECTIONS 

INDUCED LIFT CORRECTIONS 


In general, the equation for the induced lift correction at a tip 
or trailing edge due to one cons tan tr-load sector may he expressed aB 


(AL) a± = 2pV 


r ta 21 r yi(ta) r _i / , v a(x, 

l l r ^ (a ’ ta) w. 


a ii J ai 


zL 

ta)‘ 


28 / ( a— 6 ) ( a+l) ( l+t a ) d(x,y) 

1+5 J a 2 ( &— t a ) d(y,t a ) 


dy dt £ 


XB1) 


i = 1, 2, 3, 4 

y =i° 1 ^ 2 

7 i \l i = 2 


J 


The Jacobian has the value 


d(x,y) _ y-yaj 
d(y,t a ) - t a 2 


(B2) 


For i=l,2 the corrections apply at the right- and left-span tips, 
respectively, and for i=3,4 the corrections apply at the right- and 
left-span trailing edge. The limits yi(t a ) are the values of y at 
the intersections of the edges of the wing with the ray t a /p (figs. 4, 
5, and 8). For i=2 (left— span tip) it will be noted that there are 
two integrals to be evaluated in equation (Bl). This added correction 
is required at any tip which becomes a leading edge, as does the left- 
span leading edge . 

For purposes of integration it is convenient to write equation (Bl) 
in the" form 1 
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(AL) 


a i 


= 2 pv / / - 


6 ali 


cos- 1 Tii(a,t a ) - 


26 7i /(a— 6) (a+1) (l+t a ) 


(1+5) 


/ 


a 2 (8-t a ) 


£>(x,t) 

d| dt n 

d(S,t a ) a 




/(B3) 


where 


1=1, 2,3> 


6 = y-y ai 


tida) = 


yi(ia)-y ai 


(b4) 


y ai = as [p(a+6) ] — 1 «/p 2 +5 2 
ya 2 = as 1 3 ( a— B ) ] 1 «/p 2 +6 2 

y a 3 = as ^ 2+&S 

y &4 - & s (m 2 t -a 2 t ) [3(m 2t -a) (a 2t -8)] 1 J$ z +h z 


) (B5) 


and 


duA(a) duAj(a) 

^1 = ^2 = ~ da^ Ua 3 = ^4 = “!S— 


(B6) 


The values for qi and can he obtained from equations (12) and (34), 
respectively, while t) 2 and tj 4 may he obtained from ■qi and t) 4 by 
equation (Al) . 

The limits of integration for equation (B3) are 

li(la) = s t a 0i(a) [3(t a -mi^)] */3 2+ 5 2 

. Mt a ) = s t a 0 2 (a) 2 [3(t a -m 2t )] 1 ^/F+5 2 
S 3 (ta) = s t a 0 3 (a) [p(t a +5)3 A /p 2 +S 2 

i 4 (t a ) = s t a 0 4 (a) [3(t a — 6)] <^ 2 +S 2 


> (BT) 
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and the integration with respect to t a are 


tan ® ta-ip - 5 tai a - t aiA - 1 


" a 12 


13 ° a 14 


ba 23 - m l t . ta 24 “ m H ta 21 ~ ta 22 - “ 1 


(B8) 


where the functions 0-j_(a) are defined as 

(mi t +5)(a-ai t ) 


0i(a) = 


0 2 (a) - 


0 3 (a) = 


0 4 ( a ) = 


(a+6) (ai t +8) 

(m2j.— 8) (a-a 2 t) 

( a-1 -Q) (a 2 ^.— 8) 

(mi^.+8) (ai^-a) 
(mi^-a) (ai£+8) 

(m 2 ^-— 8) (a 2 ^.— a) 
(m 2t -a) (a 2t ~ 8) 




> (B9) 


The integration of equation (B3) yields for the induced lift due 
constant-load sector 


J 

to one 


(AL) ai = pVu ai S 2 Gi(a ) 


P 


} (BIO) 


i=l,2,3,4 


where the functions Gj_(a) are defined as 


Gi(a) = 0i(a) Gi(a) 


(Bll) 


The values of 0^(a) are given in equation (B9) while g^(a) are 
defined hy the following relations: 
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Si 


(a) _ 1 1 | 1 /ja 

mi t +8 m^-a mi t -a J ^ mi t 


+8) (l+a) 
+ 5 ) (mi^+l) 


82(a) = 


n^— B 


L m 2^~ a a(m2 t — 5)J 


/ (a- 5 ) (l+a) 

J (m 2t -6)(m 2t +l) 


Sc. ( a ) = 


— — + — / 

mi^+S a+B a+8 / I 


(m lt -a) (l-a) 
(l- 5 )(m lt + 8 ). 


\ (B12) 


84(a) = 


1 1 
+ 


(m 2t -a)(l-a) 


— 6 a— 8 a+B / (l—8) (m2 t -S) 


J 


For the total Induced lift (equations (16) and ( 23 )) the derivatives of 
G^(a) i=l,2 are also required. These are obtained by setting 


Gi' (a) = 20i(a)0i’ (a)gj_(a) + gi'(a)0i(a)‘ 
where the derivatives ©^'(a) are 

mit+5. 


(B13) 


01 * (a) 


(a+8) 2 


y (bi4) 


and the derivatives g^’ (a) are 

1 f (l+a) (nii|.+8)+( a +8) (mit+l) 


83.' (a) = 


8 2 ’ (a) = 


(mi t -a) 2 1 

L 2 *J (nii^+B) (l+mi^.) (a+8) (l+a) J 

1 | 

( l+a ) ( m2^.— 8 ) + ( a— 8 ) ( m2^.+l ) 

(ni2 t -a) 2 

_ 2 8) ( l+m 2 ^.) (a— 8) (l+a) 

8 

(a— 8)— 8(a+l) 


* 


) (B15) 


2a 2 (m2^.— 8) */ (n^.^.*- 8) (l+n^^) (a— 8) (1+a) 


J 
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INDUCED ROLLING-MOMENT CORRECTIONS 


The Induced rolling-moment correction for a' single elemental sector 
is defined as the moment of the induced lift of one element about the 
axis of symmetry of the wing. The general equation for the rolling 
moment for one elemental sector is, therefore. 



As in the case of the induced lift, the terms for 1=1,2 are the 
induced moments for elemental sectors at the right- and left-wing tips, 
while the terms for i=3,A are the induced moments for the right and 
left trailing edges. The values of d(x,y)/d(y,t a ) and ua*, are given 
in equations (B2) and (B6), respectively. 

For purposes of integration, equation (Bl6) may be put in more 
convenient form; namely. 



where £ and i are given in equation (b 4) while fi(t a ) and f 2 (t a ) 
are defined as 

*i(t a ) = ~tT fs(ta) = ^eT (Bl8) 
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The limits of integration are given in equations (B7) and (b 8) while the 
values for y a ^ are given in equation (B5) . The integration of equa- 
tion (B17) yields for the induced rolling moment of one element at a tip 
or trailing edge 


(AL’) ai P^ a; £ s 3 — ^2 H^(a) 


►(B19) 


i=l,2,3A 


Here the functions Hi (a) are defined as 


Hi (a) = 0i(a) hi(a) 

i=l,2,3,4 


(B20) 


where the functions Q^a) are defined in equation (B9), while the 
functions hi (a) are given as follows: 


ki(a) = gi(a) 


1 + 


Ma) f a+5 1 1-6 


m 1+ — a 2 l+m a 


:)]- 


81(a) (1—5) (a+5) 

6(l+mi^.) (mi^.+6) (mi^.— a) 


112(a) = g 2 (a) 


1 - ^ 7^ ] + ^ 7^1 


Qg ( a ) ( 1+5 ) ( a— 5 ) 

6(l+m2 t )(m2 t -5) (m 2 ^.— a) 


> (B21) 


ggi(a) = 


113(a) - 


M») = 


W 


— 6 — a. m 2 

+ (mit -1 ait) ( a+ 5) 

. (mi t -a)(a lt +6) J 

(mg-tf^t) (a— 5) 
(m2^.— a) (a 2 ^.— 6) _ 


(a— 6) ( 1+a) 
(m 2t -6) (l+m 2t ) 

83(a) 

84(a) 


J 
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Tlie values for g^(a) in equations (B21) are given in equation (B12). 
In the integration of equations (18) and (25) By parts, the derivatives 
of Hi(a) and H 2 (a) will he needed. They take the form 


(a) = 20^)0^ (a)hi(a) + 0 i (a) 2 3 h i ' (a) (B22) 


The derivatives for ©i(a) are given in equation (Bl4), while those of 
hi (a) are as follows: 


hi' (a) = gl » (a) - - - + 

6 (l+mi^.) (mi^.+5) (mi^.— a) 


0i(a) »i t +B 


3 (m lt -a) : 


8i( a ) - J 


1-6 


2 (l+m lt )(m lt +6) J 


1 r / \ n t \i * f a+6 1 1—6 

- tgi(a) 0i(a)J ~ 

3 . \mi t -a 2 l+m lt 


h 2 '(a) = g 2 '(a) - ± (H-5)(a-6)0 g '(a) + 

6 (l + m 2 t )(m 2t -6)(m 2t -^) 


> 


> (B23) 


02 ( a ) 212^—5 

3 (m^— a) 2 


1 ( 1 + 6 ) 
821 - — 


2 ( l+m 2 ^. ) (m 2 ^.— 8 ) J 


i [g2i( a ) 82(a)] ’ (-Z 2 -) - 

3 \m 2t — a/ 


I '«■« *<•»* (M;) 


62 1 * = 


(l+a) (m 2 j— 5)+(a^-6) (m 2 j+l) 
(m 2t -a) 2 L 2 J (m 2t -8) (l+m 2t ) (a-5) (l+a) 


- 1 


J 
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The values for g^a) and g^ (a) are given in equations (B12) and (B15), 
respectively. 

' INDUCED LIFT CORRECTION IN WAKE DUE TO THE 
SINGLE CONSTANT-LOAD ELEMENT 


The total lift correction in the wake due to the single constant- 
load element is given as 

i-Sk t [4 . 

i=l 1 


r 1 r 6i2 

/ / f(«pa) 


d(x,g) 

d(s >%) 


d| dt5 


where 


*-*8 x— x 5 


and the Jacobian has the value 


(B24) 


(B25) 


d(x,g) _ _PI_ 

d(l,t 5 ) t 5 2 


(B26) 


The limits t$ i refer to the rays from the apex of the trailing edge 
through the leading-edge tip. For the right span tg^ has the value 


smi ( p 2 +6 2 ) +Co&P (mx+6 ) 
s ( P 2 +5 2 )— c 0 3 ( m i+5 ) 

o. 

while the limits for the integration with respect to £ are 

^ . c 0 t 5 (m!+5) 

(tg-ffli) 8/p 2 +5 2 

= ■ 

' 12 P(t5+5) " " • 


(B27a) 


(B27b) 
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The limits tg 2 , | 21 , and g 22 for the left span may be obtained by 
applying equation (Al) to equation (B27). 


Performing the integration indicated in equation (B2k) gives 


AL 

qa 


gpV^S s2(p g +5 g ) ^ 

qa ' p l 

i=l 


Roi “ Rii 


+ 


c o g P g 

S2(p S +6 2 ) 


■? R2i 


(B 28 ) 


vhere the values of R 0 i, Rii, and R 2 i are as follows: 


where 


Roi - 


i+nit 

(l— 6) (mi^.+6) 


n Ao(q>o) / (H-mi t )(l-HBg t ) 

2 K(k) J (l-6 2 )(m H +&)(m 2t -6) 


(B29) 


A 0 (cp 0 ) - | [F(cp 0 ,k*)E(k) - F(cp 0 ,k')K(k) + E(cp 0 ,k’ )K(k) ] (B30) 


<Po 


sin -1 



(m 2t -6)(l4-mi t ) 

(l-8)(mi t +ni2 t ) 


/ ( 1-ffilt ) ( l-ffl2t ) 
(l+m lt )(l+m 2t ) 


/ (B31) 


k* = 


2(m lt +m 2t ) 
(1+mi^.) (l+m 2 ^.) 


J 


\ 

The values of Aq for various values of cp 0 are tabulated in refer- 
ence 11. The value for R 02 may be obtained from equation (B29) by 
making the transformation of parameters given in equation (Al) 


R11 = 


F(q>i,k) f 1+tpj^ 


(l+mi^.) (l+m 2 ^.) 


K(k) \(l-8)(t 5i +5) K(k)y (l-8 2 )(m lt +5)(m 2t -S) 

f ^(to) - } 


(B32) 


X 
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R 21 = (mi+8) 2 


F.(<Pi,k) 

K(k) 


I 1+tSi 
| (l+m 1 )(t 5 1 -^n 1 ) 


(l +m lt) (-*- +nl 2t) 


K(k) J (l-m 1 2 )(m lt -^n 1 )(m 2t +in 1 ) L 2 


“ Ao(<Ps) " J 2 


K(k) 

F(cpi,k). 


(B33) 


where 


Jj. = tan -1 


w 

2 £ (J 1 sin njij slnh n\>i 

n=l _ 

00 2 

1+2 Z up cos n^ cosh nv^ 
n=l 


y (B 34 ) 


1 = 1,2 


J 



(B35) 


Hi = H2 

v 1 = Jt 

V 2 = Jt 


F(cpi,k) 
' * K(k) 

^(^o^k* ) 

K(k) 

F(cp 2 ,k») 

K(k) 


(B36) 


' ) (B37) 


<Pi = sin 1 


(l+m lt )(l-t6 1 ) 

(l-^ni)(m lt +m 2t ) 


(B38) 


<P 2 = sin 1 


(m 2 ^.+mi ) ( 1+mx^. ) 
(l+mi)(mi t +m 2t ) 


(B39) 
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The values for u> are tabulated In reference 12, while k, k' , 9o> and 
tg x are given in equations (B31) and (B27). It should be noted that 
when <Pi equals Jt/2, Ji is zero. The values for R i2 and R 22 may 
be obtained from equations (B 32 ) and (B33)> respectively, if equation 
(Al) is applied. 


INDUCED ROLLING-MOMENT CORRECTION AT ROOT CHORD 


The induced rolling-moment correction corresponding to the lift 
correction in equation (B28) is 


AL' 

qa 


2 u s s3 (P g +5 g ) 

3 P qaP 




Roi ~ Rii 


+ 


i=l 


c 0 3 [mj-(-l)^S] 

T5- (p2+62) 3 


P 3 R 2 i 


C 0 3 Cmi-C-l) 1 ^] 2 3 

s 3 (p 2 + 6 2 ) 3 . P 



The expression R 31 is 


P - (“1+8 ) 2 

F(cpiA) f 

f (l+m x ) 2 

(1+mi.) 

"1 1 _ If 

(l-“i) 2 

n 3l ' 0 

(WO 2 

K(k) i 


(1-nil) (m x ^-ffli) _ 

2 


T Jf 

i A ° 


(cp 2 ) - 


K(k) 

F(cpi,k) 


x 3 - 


2 mi + 


(l-^) (Tsr-k 2 ) J 


X 


1 


(1“T 2 )(t ST* 2 ) (l“T 2 )(T 2 -k 2 )F(cp 1 ,k) L 


E((pi,k)' - 


T 2 T 1 V(l- j ri 2 )(l-k g T 1 2 T 
1 -t 2 t x 2 



(b4i) 


where 
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APPENDIX C 

FORMULAS FOR THE SECONDARY INDUCED LIFT AND 
ROLLING-MOMENT CORRECTIONS 


The corrections discussed here are those which are made at the 
leading and tip edges to correct for the residual lift and rolling 
moment introduced by the single constant-load correction in the wake. 

The formula for the streamwlse velocity for a given elemental 
sector at the leading edge of the wing is of the same form as equa- 
tion ( 19 ) , and for the right span is . 

CQS 1 (t s -m 1 )(H-t b )+(t b -m 1 )(l+t fi ) _ 2m 1 */ ( t&-mx ) ( 1+tg ) /l+tfr ’ 

(•ttr-ts) (l-Hni) ts(l+mi) // mi-tb _ 

(Cl) 

where, if x^, y^ is the apex of the element, 

(C2) 



and 


H = 3 


x b~ x 6 




dt 5 


(C3) 


(C4) 


The term in brackets to be differentiated is u(tg) of equation (24). 

The equivalent left-span solution may be obtained from equation (Cl) by 
applying the transformation of equation (Al) . In order to apply equar- 
tion (Cl) at the right tip, it is necessary to replace mi by —5 ana 
at the left tip by 5. It will be noted that at the right tip the 
second term of equation (Cl) is imaginary, while at the left tip it is 
real. Therefore, it is necessary to retain the second term in the square 
brackets for evaluating the correction along the left span. 
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Induced— Lift Corrections 


The secondary correction to the lift for a single element at the 
leading edge of the right span takes the form 


— (AL) dtp =~^{ r 1 f f cos 1 »l( t, b» m i) “ 

dtg L.E. n x J Q 

^ , 2 ^ dtt> + 

% 

Ujj [cos -1 T)(t-b,mi) - 

•T.i(tb,mx)] ~ dft dtb j (C5) 



vhere In and I 21 correspond to the intersection of the rays t"b with 
the trailing and tip edges, respectively, of the wing plan form, and 
tb! is the ray tb passing through the trailing— edge tip. The expres- 
sion in brackets denotes the bracketed expression in equation (Cl) . The 
limits of integration have the values 


| - c o t bfo 1 - ,f8 )( t 9r ni t) 

11 */p 2 +5 2 (m lt — tb)(t 5 -mi) 

[s(fl 2 +5 2 ) (tg-mi')— Pc 0 (t5+5) (mi+5)] 

0 (tb+8) (tg-mi) «/p 2 +5 2 

3 c o‘ fc 6( m i + ^) (mx t +S)-sm lt (3 2 +S 2 ) (tg-mx) 
Pc 0 ( mx+8 ) ( mx t +5 )-s ( |3 2 +S 2 ) ( ts-mx ) 


(C6) 

(C7) 


(08) 


while the left-span correction may be obtained from equation (C5) by 
applying equation (Al) . 

The secondary correction to the lift along the right-span tip for a 
single element is 
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1 tip dt5 = f /' “>> Cc ° S 1 “ 


dt 5 


(AL) 


-i w o 


T(t b ,-6)] -L d| dt b 
tb 


(C9) 


The limit of integration | 31 refers to the intersection of the ray t b 
vith the trailing edge, and has the value' 


st b (t 6 -^ni t ) ,/p g +5 g 
0 ( 1 5+8 ) ( mi b — tb ) 


(CIO) 


The lefi>-s pan-tip correction may be obtained by applying equation (Al) 
to equation (C9) . 

The total secondary correction to the lift is obtained by summing 
the lift due to all the elemental sectors along the leading and. tip 
edges. The correction along the tip includes an initial sector of con- 
stant magnitude F[cp(t5 1 ) ,k/k(k)]u 5 . It should be pointed out that due 
to the different orientation of the elemental sectors along the tip and 
leading edges there is an overlapping of sectors in the region bounded 
by the extensions of the leading edge and the ray tg^ from the apex of 
the trailing edge through the tip of the leading edge. An approximate 
method of correction for this overlapping is to let the initial tip 
sector be bounded by the extension of the leading edge rather than by 
tg^. Then the total secondary correction to the lift for the right span 
is obtained as 


(**)i . i { / 




L 


mit (AL) dt 6 + 
t 5l dtg tip 


2pV3 u F[cp(tg 1 ),k] p 6 P^3ir_ x (mi+5) (l+tb)+(tb+S) (14-mi) 

it 5 K(k) J 0 L '( t b — fflj_ ) ( 1—5 ) 

25 V(m^5)(lt m3j. j 1 * d | dt t \ (Cll) 

mi (1—6) J -(5+tg) J ^ J . 

vhere 

✓ 


t 5i 


mis(P a +-6 2 ) + c n 5ft(mi+5) 
s^ a +5 a ) — Co3(nii+6) 


(C12) 
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The total secondary correction to the lift at the left span may he 
obtained from equation (Cll) by transformation (A5). If the Mach cone 
from the apex does not cross the leading edge, then it is necessary to 
evaluate only the second integral of equation (Cll), since the limit 
t 5i becomes unity. 


Induced Rolling Moment 


The secondary correction to the rolling moment for an elemental 
sector at the leading edge is 




/'■“■* '•-■■'--.I- 


— 1 o 


TiCtb,*!)] | d | dtb + 


n m 1 . n^zx 

J J % f cos - ^ Ti(t t ,m!) - ^(t^mx)] x 
tb x o 

fS fxi(th)+(S Si )f2i(tt>)3 
*b*. 


i dtb 


} 


(C13) 


where fn(tb) and fai(tb) are defined as 


f 11 = 


f 2 i - 


tb+8 

tb 


t 5 +5 

t 5 


(C14) 


and 


5 81 = 


c 0 t6(m!+5) 
•J p 2 +5 2 'tg-ni! 


(C15) 


The limits of integration are as given for the lift in equations (C6), 
(C7), and (c8) . 
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The sec ondar y correction to the rolling moment for a single element 
along the right tip takes the form 


(AL' ) , dt 8 = g P v P. 2 — r r [cos-1 T|(tb,-8)- Tft^-b) ] x 
dt 5 1 tip J J 


[|fn(tb)~ > ~(^Si)f2i(th) 1 


t d| atb 


The value for in equation (Cl6) is 


sts Vf3 2 +5 2 

l& i = p(t 8 +5) 


(Cl 6) 


(C17) 


■The left— span correction is obtained from equation (Cl6) by transforma- 
tion equation (Al) . The limit of integration Ig-^ is the same as pre- 
viously given in equation (CIO). 

The total secondary rolling-moment correction for the right span is 


AL* 

qa 


qa 



*81 _d 

dt 8 


& - dt & + 


, nil 


L.E. 


77- (AL’ ) dt 5 + 

dt B Hip 


2pVP 2 F[cp(t& 1 ),k] 

u 6 


X 8 X 

-~8 -,£31 


I f 

J — 1 J o 


*y^ 2 u K(k) 

[|fii(tb) + (I 8i)^2i(tb)] -^2 dl dtb} 


C(mi,S,tb) X 


(Cl8) 


where the limit t Bl is as given for the lift in equation (C12) , the 
value for S 8l is given in equation (C17) , and C(mi,S,tb) is the 
bracketed expression in the last integral in equation (Cll) . The left- 
span correction is again found by applying equation (Al). If the Mach 
cone from the apex of the trailing edge does not cross the leading edge 
then all terms of equation (Cl8) but the second vanish. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Oct. 30, 1950. 
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TABLE I.- CALCULATED ROLLING MOMENT FOR WINGS 1 AND 2. 


Source of moment 

L , Wing 1 

— per degree of 
yaw 

t Wing 2 

— per degree of 
yaw 

Basic 

10.70 ft® 

73. 4# 

2.33 ft® 

97.00 

Tip effect 

.04 ft 3 

•30 

.07 ft 3 

• 3.00 

x Wake from apex 
of trailing edge 

4.32 ft 3 

29.6 0 

— 

— 

Oblique wake 

— 

— 

— 

— 

^■Secondary 

-.48 ft 3 

-3.# 

— 

. 

. Totals 

14.58 ft 3 

1000 

2.40 ft 3 

1000 


qa 
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Figure 3~ Velocity distribution u^(o) for a yawed triangular wing 


X 


Section f-f 


Figure 4. 



— Typical sector (shaded) and the boundary conditions 
for the right- span -tip correction. 
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Figure 7.— Typical sector (shaded) and the boundary conditions 
for the correction from the apex of the trailing edge . 







X»|g 



Figure 8.— Typical sector (shaded) and the boundary conditions 
for the right-span oblique trailing-edge correction. 
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of attack with yaw angle ^ for M=V2\ 












